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North Carolina State University 

Let {sij — 1,2, .. .} consist of i.i.d. random variables in C with 
Esii = 0, E|siip — 1. For each positive integer A^, let — Sk(N) — 
{sik,S2k,-- ■,SNk)'^, l<k<K, with K = K[N) and K/N -> c> as 
A'^ — » oo. Assume for fixed positive integer L, for each N and k < K, 
OLk = (afc(l), . . . ,OLk{L)Y' is random, independent of the Sij, and the 
empirical distribution of (ai, . . . , oik), with probability one converg- 
ing weakly to a probability distribution H on C^. Let (if. — f3i.{N) = 
{a^Wsl, ak{L)slf and set C = C(iV) = (l/N) J2k=2 Pkl^l- Let 
cr^ > be arbitrary. Then define SIRi = [1/ N)Pl{C + cr'^ I)'^ P^, 
which represents the best signal-to-interference ratio for user 1 with 
respect to the other K — 1 users in a direct-sequence code-division 
multiple-access system in wireless communications. In this paper it is 
proven that, with probability 1, SIRi tends, as ^ oo, to the limit 

(^i-'i^oii{()o.i{i')ai^ti , where A — {aiji) is nonrandom, Hermitian 
positive definite, and is the unique matrix of such type satisfying 
A = (cE ""' -f a^lL)~^, where a € has distribution H. The 
result generalizes those previously derived under more restricted as- 
sumptions. 

1. Introduction. This paper examines the mathematical properties of a 
quantity fundamental in analyzing the performance of a particular scheme 
used in wireless communications. The scheme, known as direct- sequence 
code-division multiple- access (or DS-CDMA), currently in use, effectively 
handles many users by taking into account the manner interference inter- 
acts when a particular user's information is being decoded. It is achieved by 
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assigning to eacli user a vector of high dimension, cahed a signature sequence. 
Suppose there are K users and L receive antennas. Let be the dimension 
of the signature sequences, and denote by e the signature sequence 
assigned to user k. At a particular instant of time let X/. G M denote the 
value transmitted by user k having transmit power G M"*", and let 7fc(^) 
denote the fading channel gain from user k to antenna ^. It is assumed that 
the XkS are independent standardized random variables. With W(£) G 
denoting noise associated with transmission to antenna i, entries Wi{i) i.i.d. 
across i and i, mean zero and E|Wi(^)p = u^, the data recorded at antenna 
£ is modeled by 

K 

k=l 

Letting Y = [Y(l)'^, . . . , Y(L)-^]^ G C^^, the goal is to capture the transmit- 
ted Xk for each user in a linear fashion, that is, by taking the inner product 
of Y with an appropriate vector G , called the linear receiver for user 
k. For user 1, Xi = c^^ Y is the estimate of transmitted Xi. 
The output signal-to-interference ratio 

|ctSl|2 



f^'l|ciP + Ef=2|c|Sfc|2 

associated with user 1 is typically used as a measure for evaluating the 
performance of the linear receiver. Here 

sk = Vnhk{i>L...,7k{L)sif. 

It turns out that the choice of ci which minimizes E(X — X)'^ (the minimum 
mean-square error) also maximizes user I's signal-to-interference ratio, the 
latter taking the value 



\k=2 ) 



Si, 



where / is the NL x NL identity matrix. It is this quantity which is the 
focus of this paper. 

In [4] properties of SIRi and their dependency on the 7fc(£)'s, T^'s, o"^, 
L, and when the latter two values are large, are explored by proving 
limiting results, as and K approach infinity with their ratio approaching 
a positive constant, under the assumption that the s^'s are randomly gener- 
ated (which is usually done in practice). They are independent i.i.d. random 
vectors containing i.i.d. mean zero entries, independent of the 7fc(£)'s and 
Tfc's. The results allow for analysis of performance in various situations de- 
pending on the location of the users with respect to each other and the 
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antennas. Additional assumptions are imposed on the 7fc(£)'s and the T^'s. 
Throughout [4] it is assumed the jk {(■) 's are independent and circularly sym- 
metric (i.e., the argument of each 7fc(£) is uniformly distributed on [0, 27r)), 
and the entries of each s^, are mean zero complex Gaussian with variance 
l/N [i.e., they are of the form Zi + iZ2 with Zi, Z-z i.i.d. iV(0, l/(2iV))]. The 
Tfc's are allowed only to depend on |7j(^)| for all k,j,i. Two scenarios de- 
pending on the location of the antennas are considered. One scenario places 
all the antennas near each other, the other allowing them to be located 
anywhere. Theorem 1 in [4] applies to the former case, the proof of which 
requires, for each k, 7^(1), . . . ,^k{L) to be identically distributed. Also, T^, 
as a function of the |7j(£)|'s, is assumed to be symmetric with respect to 
the antennas, in the following sense. For any permutation vr on {1, . . . ,L}, 
we have for each k, 

n{T{n{l)), T{7t{L))) = Tfc(r(l), . . . , r(L)), 

where T{£) = (7i(^), . . . , 7fc(^)). These two assumptions are lifted in Theo- 
rem 3, which would not be realistic when antennas are not placed in one 
location. However, an additional assumption is made, namely, that there are 
for each user independent signature sequences going to the L antennas, that 
is, Sfc takes the form 

with s^j^\ . . . , s^f^^ i.i.d. As pointed out in [4], this "completely random se- 
quence model is not physically realizable." 

The purpose of this paper is to prove limiting results on SIRi substantially 
more general than those found in [4]. The main result is the following: 

Theorem 1.1. Let {sij = 1,2, . . .} be a doubly infinite array of i.i.d. 
complex random variables with Esu = 0, E|siip = 1 (we will from henceforth 
call standardized). Define for k = 1,2, . . . , K si^ = Sk{N) = (sifc, •S2fc) • • • > SNk)'^ 
We assume K = K{N) and K/N ^c>0 as N ^ oo. For each N , let 
jf.{i) = ^^(^) g C, Tk = T^ £R^, k = l,...,K, i = l,...,L, be random 
variables, independent of si, ... , sk's. Let for each N and k, 

ak = ak =VT'kilk{l),---,7kiL)f. 

Assume, almost surely, the empirical distribution of cx\, . . . ,C)lk weakly con- 
verges to a probability distribution H in C^. 

Let (3k = (3k{N) = Vn{jkil)sl, . . .,7k{L)slf and 
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Define 



then, with probability one, 

L 

lim SIRi=Ti 

W— >oo 



where the Lx L matrix A = (a^ £/) is nonrandom, Hermitian positive definite, 
and is the unique Hermitian positive definite matrix satisfying 



where cx G has distribution H and II is the L x L identity matrix. 

Clearly SIRi defined in this theorem is the same as the one initially in- 
troduced, the only difference in notation being the removal of the scaling by 
in the definition of the s^'s. 

Let ai denote the ith entry of the random vector cx having distribution 
H. Under the independence and circularly symmetric assumption on the 
7fc(^)'s and the independence of their angles and the T^'s, it follows that for 
£^ i' and positive ai, . . . , a^, 

(1-2) E ^ ' =0. 

With just this additional condition we have the following: 

Corollary 1.1. Under the conditions in Theorem 1.1 and (1.2), the 
limiting A = diag(ai, . . . ,0^), where the a/s are positive satisfying 

(1.3) a ^ 



' cE(|a,|2/(l + Eia,|a,|2)) + c72- 

Corollary 1.2. Suppose the conditions in Theorem 1.1 are met except, 
for the limiting behavior of the cx-k s, it is only known that: 

1. the empirical distribution of 

(1.4) Tfc(|7fc(l)|2,...,|7fc(L)|2)^, 2<k<K, 

converges almost surely in distribution to a probability distribution G in 
M^, and 
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2. for £' and positive oi, . . . 



,aL 



1 



K 



(1.5) 



K-1 



E 







k=2 



almost surely, as N ^ oo. 

Let ((5i, . . . , 6l)^ E denote a random vector having distribution G. Then 
the conclusions of Theorem 1.1 and Corollary 1.1 hold, with each ja^p in 
(1.3) replaced by 6i. 

Theorem 1.1 frees up conditions on the 7fc(^)'s, and overall dependence 
between them, the T^'s, and the antennas. Moreover, the result allows for 
more general (and realistic) assumptions on the generation of the s^'s, per- 
mitting their entries, for example, to be just ±1, which is typically done in 
practice. 

Thus, under the general assumptions in Theorem 1.1, various scenarios 
can be analyzed and compared. In applications the empirical distribution of 
the Qfc's would typically be used for H, the matrix A thereby satisfying 



Although there appears to be no explicit solution to (1.1), it will be shown 
that A can be computed numerically by iteration of the right-hand side of 
(1.1), provided the eigenvalues of the initial choice in the iteration lie in a 
certain closed interval in (0, oo). 

The conclusion of Corollary 1.1 is the same as that of Theorem 3 in [4], 
where only the a.s. convergence in distribution of the empirical distribution 
of (1.4) is assumed. But recall Theorem 3 also assumes for each user different 
signature sequences for each antenna. The extent of the results in [4] are 
confined to diagonal ^'s due to the assumptions imposed [essentially the 
spherically symmetric assumption on the 7fc(^)'s], clearly a special case of 
the general conditions assumed in Theorem 1.1. 

It is remarked here that the assumption of the Sij coming from a doubly 
infinite array can be replaced with Sij = Sij{N), 1 <i < N , 1 < j < K , with 
no dependency assumptions for different N, provided E|sii|^ < oo. Indeed, 
it will be seen in the beginning of the proof of Theorem 1.1 that the dou- 
ble array and finite second moment assumption is needed only when the 
strong law of large numbers is invoked on sums involving | Sij \ ^ , the alterna- 
tive assumptions yielding the same conclusions with the aid of Lemma 2.10 
below. 

Theorem 1.1 only provides limiting properties of the signal-to- interference 
ratio with respect to one user. The last section of this paper will address the 
issue of uniform convergence of all the K SIR's. 
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The proofs of these results will be given in Sections 3-7, with basic math- 
ematical results needed in the proofs presented in Section 2. 

Note. After submitting this paper, the authors came upon a result sim- 
ilar to Theorem 1.1, announced in a conference paper, without proof [3]. 
In that paper it is claimed that a proof is given in another paper, sub- 
mitted for publication. In [3] the Sjj's need not be identically distributed, 
nor come from one doubly infinite array of variables, but it is assumed 
E|sjj|^ = E\sij{N)[^ < N'^~'^ for some 7 > 1. Moreover, the limiting distribu- 
tion H is assumed to have bounded support. The conclusion has convergence 
in mean square. 

2. Basic tools. This section contains properties of matrices, a classic 
fixed point theorem, and some probabilistic results, needed in the proof of the 
above statements. Throughout, / will denote the NL x NL identity matrix. 
For arbitrary dimension n, In will denote the n x n identity matrix. For 
any rectangular matrix X, vecX will denote the column vector consisting 
of stacking the columns of X on top of each other, first column on top, last 
on bottom. Spectral norm on matrices and Euclidean norm on vectors will 
be denoted by || • ||. 

Lemma 2.1. Let cj^ > 0, B, A n x n matrices with B Hermitian non- 
negative definite, and x € C" . Then 



Proof. The identity follows from (Z? + xx*)"^x = D~^x ^^^,^_i^ , true 
whenever nxn D and D + xx* are both invertible. Write B = J2 ^i^i^i^ its 
spectral decomposition. Then 



tr((S + XX* + -{B + a^I)-^)A\ 

_ x*{B + a'^I)-^A{B + (j^/)"^x ^ 
l + x*(5 + cr2/)-ix - 



< 



a- 



A 



■2 



l + x*(S-hcj2)-ix 



^11 ||(5 + Cj2/)-lx 



- l+X*(5-f-cr2/)-lx 



E(l/(A. + a2)2)|e*x| 



l + E(l/(Ai + a2))|e:x| 



A 




□ 



The next lemma is easily verifiable. 



SIR OF CDMA SYSTEMS 7 

Lemma 2.2. For any matrix A N x K and a"^ > 0, 

{AA* + cj2/^)-i = - A{A*A + a^/^)"^^*). 

The following lemma is a direct consequence of the previous identity. 

Lemma 2.3. Suppose Ai,. . . ,Al are N x K, and cj^ > 0. Define the 1,1' 
block of the NL x NL matrix A by Ai^£i = AiA'^, and, splitting (A + a'^I)~^ 
into L? N X N matrices, let {A + a'^I)j^, denote its i,£' block. Then 

{A + a^I)j} = (^6i,i,I^ - A^ A^Ae + a^lK^ A}}j . 
Lemma 2.4. Given Ai,...,Al are N x K and zi,...,Zi € C with 



Y^AzA ^E^^^^ 



where represents the partial ordering on Hermitian nonnegative definite 
matrices. 

Proof. For any b G C^, we have by two applications of Cauchy-Schwarz, 



b* ( ^ezi I ( E A}zi ]h = Y, h*AiA}hziZi 



< J2{^*AeA}hy/\h*Ae>Al,h)'/^\ze\ \ze\ 



= (Y(^*A,A}h)'/'\z,\j 

< E l^^l' E bM,^|b = b* (y AeA}) b. 



This proves the result. □ 



Lemma 2.5. For Ai, . . . , A^, A, cr^ in Lemma 2.3, the L x L matrix 
{^{A + a'^I)jj,) is positive definite with smallest eigenvalue bounded below 
by 

trf^A^^I + a^/^^ . 
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Proof. For zi,.. .,ZL&C,J2e = 1' we have by Lemmas 2.2-2.4 the 
smallest eigenvalue of (tr(A + a'^I)jj,) is bounded below by 

^ tr(A + a'^l)jj,zeze' 

= 0-"^ tr Se^e'ZezplN - X! X! "^^"^^ + ^^^^ M'ZiZi>i \ 
= tr {in - Atz^ A} A, + a^Ii^ Alz^j j 

X Alz^ A,z^ + a^Ij^ Alz^ ^ 



tr [Y,AizA[Y,A}zA +cj'^In] > tr + a^/^ 



-1 

■ □ 



For A = (aij) m x n and B p x q, the Kronecker product of A and B, 
denoted by A ^ B, is the mp x nq matrix, expressed in blocks oi p x q 
matrices, the i,j block being aijB. We will need the following, which is 
Lemma 4.2.10 of [5]. 

Lemma 2.6. For Amxn, Bpxq, Cnxk and D qx r, we have 
{A ®B){C®D) = {AC) ® (BD). 

The following is needed to prove Corollary 1.1. 

Lemma 2.7 (Schauder fixed point theorem [7]). If A is a convex, com- 
pact subset of a Banach space X and g:A^ A is continuous, then g has a 
fixed point in A. 

The next result is one on the eigenvalues of the expected value of the 
Kronecker product of two random matrices. 

Lemma 2.8. Let A = (oij) = (ai, . . . , a„) (mx n) and B (hx g) he two 
random matrices, the entries having hounded second moments. Then 



EA®B\\ <min(i/||Eylyl*|| ||ES*B|U/||E^*^|| \\EBB*\\). 
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Proof. For any hxm X = (xi, . . . , x^) and gxn matrix Y = (yi, . . . ,y„) 
with trXX* = tryy* = 1, we have, using Cauchy-Schwarz, 

\\{vecX)*[EA(^B]vecY\\^ 



^^Ea,jX*Byj 

i=ij=i 



E EaJX*i?y, 



< EE[ll^a,||||5y,||] < K:EV2||Xa,fEV2py,| 



\j=i 



\j=i 



< Ell^aj ||2 ^ E\\Byj f = ti{EAA*X*X) tT{EB*BYY*) 
j=i i=i 

< ||E^^*|| \\EB*B\\. 
Notice that 

n 

J2 ^aJX*Byj = tr EA^X*BY = tr EB'^XAY'^, 
i=i 

so we also have 

||(vecX)*[Eyl0 5]vecyf < \\EA*A\\ \\EBB*\\. 
The truth of the lemma follows. □ 

The next result, which is Lemma 2.7 in [2], constitutes the main contri- 
bution of randomness to Theorem 1.1. 

Lemma 2.9. For X = (Xi, . . . i.i.d. standardized entries, C nxn, 

we have for any p>2, 

E|X*CX - trC|P < Kp{{E\Xi\^trCC*)P/'^ + E|Xi|2Ptr(CC*)P/2), 

where the constant Kp does not depend on n, C , nor on the distribution of 
Xi. 

The last two results provide conditions guaranteeing the strong law of 
large numbers. 

Lemma 2.10. ([6]). ForXi,X2,... i.i.d., letSn = Xi-\ hX„. For 

t >1, the joint conditions E|Xi|* < oo and EXi = b are equivalent to the 
condition 

oo 



71=1 





n 



> e] <oo 



for every e > 0. 



10 



Z. D. BAI AND J. W. SILVERSTEIN 



Lemma 2.11 (Lemma 2 of [1]). Let {Xij, i,j = 1, 2, . . .} be a double array 
of i.i.d. random variables and let a> ^, P >0 and M > be constants. Then 
as CO, 



max 

j<Mnf^ 



n 



i=l 







a. 8. 



if and only if the following hold: 

E|Xn|(i+^)/"<oo 

and 

EXn, ifa<l, 
any number, if ol> 1. 

3. Proof of Theorem 1.1. For the remainder of this paper we write 
Tk-ikii) as au{l). Write 

SIRi = ^ E «i + ' 



where (C + cr'^I)^}, N x N,is the £,£' block of the NL x NL matrix (C + 
a'^I)"^. Some of the NL x NL matrices below will also be viewed in block 
form: L'^ N x N matrices. 

We begin by truncating and centralizing the entries of si. For each A'^, 
define 

Snl{N) = S„,l/(|s„i|<(l/3)logAr) " Es„l/(|s„i |<(i/3) log TV) , 

Si = Si (N) = {sn{N),S2i (iV), . . . , sm{N)f, 

^1 = ^i(iV) = (ai(l)sf , . . . , ai(L)sif , 

and SIRi = i^;(C + a2/)-i^i = ^ E^/' «i WailO^t (C + cj2/)-1si. 
We have 

ISIRi - SIRil < a~^N~\\\l3, - + 2\\(3,\\ \\(3, - M) 



^iV-i(||si-sif + 2||si||||si-si||). 



We have by the strong law of large numbers (SLLN) N ||si| 
iV — > oo. Also, since 



1 a.s. as 



N 



Es„i/(|s„i|<(i/3)logAr)| 



N ^||si -Slip -iV El^"lA|sni|>{l/3)logAr) 
n=l 

N 

< iV~^2 ^ \Snl\^ I{\s„i\>{l/3)logN) +2|Esii/(|j,-^^|<(i/3)logAr)| , 



n=l 
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TV 



limsupiV ^||si-sif< lim iV ^ V 2|s.„i/(|s„j>m)| 



n=l 



2E|si 



a.s. 



'Ill J(|sii|>M) 

Since this last expression can be made arbitrarily small by choosing M 
sufficiently large, we have 

|SIRi -SIRil ^0 a.s. 
Let si = si(iV) =si(7V)/(E|sii(iV)|2)i/2, n ig ^lear that E|sii(A^)|2 ^ 1 
as — > oo. Therefore, 



1 







as N ^ oo. Since the entries of si are i.i.d. standardized and, for all large 
A^, bounded by logA^, we have by Lemma 2.9, for any N x N A, 

E\N~^sl Asi - N~Hv A\^ <K4\\ A f (log NfN"'^, 

which is summable for A bounded in norm. Since we have \\{C + a^I)jl\\ < 
(7~^, we conclude that, with probability one. 



SIRi - iV-i ^ ai (£)qi {£') tr{C + (J^I)^} 



11' 



0. 



Similar to what is argued above, we have by SLLN 
^ N K 

Ji^^]^ E E \snk\^hKk\>iogN) = a.s. 

n=l k=2 

It is straightforward to verify the existence of a nonrandom sequence {oat}, 
of positive numbers increasing to infinity, satisfying 

o N K 



(3.1) 



TV) 







n=l fc=2 

almost surely, and 

(3.2) a7vEsiiI(|sj^l<iogAr) ^ 0. 

We may assume cat < log A^. 

Define /Ci = IC{N) = {k e {2, . . . ,K} : max£ \akii)\^ < a^} and /Cs = {2, . . . , 
K} — ICi. Since the empirical distribution of converges weakly a.s. to a 
probability distribution in C^, we must have, with probability one, 

#]C2 =Kn = o{N). 
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Write /C2 = {/ci, . . . , kKf,}, Cq = C, and Cj = Cj-i- (3i,.(3l. ,j = l,..., Kn- 
Let le'/ denote the NL x NL matrix consisting of the N x N identity matrix 
in the block, zeros elsewhere. Then, using Lemma 2.1 for any 

\N"^ trie + a\} - N-^ tr(C'i^^ + a^I)i}\ 

J2tr{Cj^i + a^iy^Ie,e-N-Hv{Cj+a^iy^Iei 
1 ^ N-'(3l(Cj + aH)-H,,4Cj + a^I)-'P, 



l + iV-i/3t.(C,+a2l)-i/3, 



< 



Kn 



a.s. 



Therefore, we may assume each q^^. =0, and that m.ayi2<k<K,e \ (^k{^)\'^ ^ cln- 
We truncate and centralize the entries of each Sk, 2 < k < K , in the same 
manner as performed on Si, and call it Sfc. The corresponding /3j, and C are 
denoted by and C. For any £,£', we have 

\N-\tiiC + a^I)jl - tT{C + a^I)lj,)\ 



K 



tr((C7 + a'l)-\(3,(3l - Mk){C + a^lY^h 

k=2 
1 ^ 



k=2 



PliC + a^I)-'h',e{C + a^ir'p, 



<^,T.\\(^k-Pkf+mk\mk-Pk\ 

k=2 
K 



K 



1 



5] 5] |afcW|2(||sfc - Sfcf + 2||sfc|| llsfc - Sfcl 



< 



k=2 I 
K 



K 



1/2 / K 



.k=2 



\k=2 



\k=2 



By SLLN, we have 
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and from (3.1) and (3.2), 

iV2 



2 K 

N II ~ ||2 

2^ ||Sfc -Sfcll 

k=2 



knfcl ^(|s„fe|>(l/3)log7V) 

n=lfc=2 

2 I 1 2 

+ 2a^|Esii/(|s-^^|<(i/3)iogAr)| ^0 a.s. 

Therefore, 

N-\tr{C + (t2/)-i - tr{C + a^/)"!,) ^ a.s. 

It is easy to verify 

A^"i(tr((E|sn|2)-iC' + a2/)-i - tr((7 + a2/)~i ) ^ 0. 

Therefore, returning to the original notation, we may replace the doubly 
infinite array and assume for each A^, Snk = Snki^), I < n < N , 2 < k < K , 
i.i.d. standardized random variables bounded by logA^. The quantities = 
Sfc(iV), Pf, = Pi^{N), and C = C{N) are defined accordingly. 

Define = C - (l/7V)/3fc/3^. Select e € (0, 1/10). Applying Lemma 2.9 
to each of s^, (C(fc) + a'^)^^,, 2 < k < K , I < < L, with p = 5, along with 
standard arguments using Chebyshev's and Boole's inequalities, together 
with Lemma 2.1, we have 



(3.3) max 

ke{2,...,K} 

e,e'e{i,...,L} 



— (s^(C(,) + a^I)llsk - tv{C + a^I)-! ; 



a.s. 



Define the L x L matrix B_= {b^^i) with 

1 ^ akinM^) 



he' - ]v E Y 

k=2 



+ (l/iV)/3^(C(fc)+a2/)-i/3, 



and define the NL x NL matrix B in terms of the Kronecker product: 
B = B®In. We have {B + cj^/)-! = {B + a^h)'^ In- Denote the £J' 
entry of {B + a^lL)~^ by 6^,^/. 
We write 

1 ^ 

C + - (B + a^I) = - ^ /3,/3^ - i?. 

fc=2 

Taking inverses on each side, we have 
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K 



N 



k=2 



Multiplying on the right by Ii'/, taking traces and dividing by N, we get 
iV-i tr{B + - tr(C + 

_ 1 ^{_l/N)(3UC^k)+aH)-^Ie'AB + c7^I)-'Pk 



N 



k=2 



l + (l/iV)/3^(C(fc)+a2/)-% 



- N-^ tr B(C + a^I)~^Ie',e{B + a^/)"! 



N 



E 

fc=2 



i + (i/iVM,(qfc) + a2/)-i/3, 



E 



.k=2 



E 



N 



l + (l/iV)/3I,(C(,,) + a2l)-i/3, 
tr 5^,^^[(C + a'^iyHe^B + 

l + (l/iV)/3^(C(fc)+a2/)-i/3fc 



.k=2 



tr 5^/ ^(C + a^I)jl, (B + a^/),-;, 



1 ^ 

^ iv ^ T 



bee'C^kiDakil') 
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Using (3.3), the fact that the 6^£''s are bounded by cj~^, and noticing that 

we immediately get 

(3.4) \bi^e' - N-^tr{C + a^I)j}\^0 a.s. 

Notice Ci^i' , the i, £' block of C, can be written as 

where S" = (s2, . . . , s^^) and ^(£) = diag(a2(^), . ■ • , afc(^)). Therefore, from 
Lemma 2.5, the Lx L matrix (tr(C + cj^I)^^,) is positive definite. 
Using this and (3.3), we have 



(3.5) 



, 1 ^ ak{f)ak{i) 



N 

'.,£' k=2 



tr(C + cr2/)7i| ^0 a.s. 



From Lemma 2.5 the smallest eigenvalue of = {N ^ tr(C + o"^/)^ j,) is 
bounded below by 

(3.6) l.tr(j^Sj2AmiirS*+a'lr?j . 

This quantity is the Stieltjes transform of the empirical distribution of the 
eigenvalues of 

^sY^Amiiys* 

£ 

evaluated at — cr^. We have, with probability one, the empirical distribution 
of the diagonal entries of 

converging weakly to a nonrandom probability distribution. Therefore, from 
[8], we see, with probability one, the empirical distribution of the eigenval- 
ues of (3.6) converges weakly to a nonrandom probability distribution, and 
consequently, (3.6) converges a.s. to a nonrandom positive number, say, m. 
Therefore, 

(3.7) liininf Amin^Af > a.s. 
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Consider a realization in which (3. 4), (3. 5) and (3.7) hold and the empirical 
distribution of {a2, . ■ ■ ,cxk) converges weakly, where a G denotes the 
random vector having distribution H. Let {Ni} be a subsequence in which 
each A^~-'^tr(C + a'^I)^^, converges, say, to ae/' for S {!,..., L}. Let 
A = {a£^£i) and 6 = infjv^ Ajnin(^7Vi) > 0. We have for z G C^, 



zz 



^ \\zf ^ 1 

- l + A„,in(A)||z||2 - 6- 



l + z*Az 

Therefore, by the dominated convergence theorem, along {Ni}, 



—t- 

' - 1 1 



K - 1 ^ 1 + Qt^ctfc 1 + aMa 

and since 



zz zz 



l + z*^Arz 1 + zMz 
we have by (3.5), along {N^}, 

„ _ cxa 



<hAN-A\\, 





1 + a*Aa' 

So A satisfies (1.1). The next section shows that only one Hermitian positive 
definite A will satisfy this equation. With this fact we have, with probability 
one, A]\i converges to a nonrandom Hermitian positive definite Lx L matrix 
A satisfying (1.1). 

4. Proof of uniqueness. Suppose A and A are two different L x L Her- 
mitian positive definite matrices satisfying (1.1). Then 

~ Aaa* Aa* (A - A)cx 

A-A = cE- 



(1 + a*Aa){l + a*Aa) 

Multiplying A~^^'^ on the left and A~^^'^ on the right, we obtain 

^1/2^-1/2 _ ^-1/2^1/2 ^ ^^ A^/^aa*Ay^a*{A -A)a 

(1 + a*Aa){l + a*A(x) 

„^*^*Ml/2^-l/2 _ ^-1/2^1/2)^ 

= cE , 

(1 + a*Acx){l + a*Acx) 

where rj = A^l'^a. and r) = A^l'^a.. Write /x = vec(A^/2^-i/2 _ ^-i/2^i/2)_ 
With the aid of the Kronecker product, we can write the above equation as 

(4.1) ^ = cE ^. 

{l + a*Aa){l + a*Aa) 
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Using Lemma 2.6, we have 



' (1 + a!MQ)(l + a*A(x) 



cE- 



cE 



l + a*Aa l + a*Acx, 



and, since 7^ 0, this matrix has an eigenvalue equal to 1. By Lemma 2.8, 
its largest squared eigenvalue cannot be greater than 



cE 



1 + cx*AcxJ 



cE 



We have 



and since 



cE 



1 + a*Aa, 
Ay^aa*Ay^ 



cE 



l + a*Aa, 

^1/2 

Ota 



1 + a*Aa 

is nonnegative definite, we have 



{l + a*Aay 
A^/^aa*Ay^cx*Aa 



cE 



2 ^^V2oaMV2 
-< cE 



l + a*AaJ " l + a*Aa 

= A'/\A-' - a^lL)A^/^ = Il- <y^A, 



the eigenvalues of which must all be less than one. The same result applies 
for the other matrix involving A. Therefore, the matrix in (4.1) cannot have 
an eigenvalue equal to one, a contradiction. So we conclude that there is 
only one Hermitian positive definite solution to (1.1). 

5. Convergence of iterations. Let f{A) denote the right-hand side of 
(1.1), considered as a mapping of the set of Hermitian positive definite ma- 
trices, which we will denote by 7i. Clearly / maps 7i into itself with largest 
eigenvalue not larger than a~'^. We proceed in finding a positive h < cr"^ 
such that / maps 

n[b, o-"2] ={Aen: all eigenvalues of A lie in [b, a"'^]} 

into itself. Notice from the dominated convergence theorem 

|2 



g{a) =cE 



a\\a\ 



1 + allaiP 



is continuous and nondecreasing for a G [0, o'"^] with g(0) = and g{a~'^) > 1. 
Therefore, there exists a G (0, cr"^] for which g{a) = 1. We claim a suitable b 
is a/(c + 1). Indeed, suppose the eigenvalues oi A ^Ti are contained in 
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[a/(c + I), a If a = Amin(^) > o, then using the fact that ||Ei?|| < E||i?| 
for any random matrix B, 



Amax(/-'(^))=^' + cA„ 



aa 



1 + a*Aa 



<a^ + cE 



a 



1 + allalp 



< 



1+c 1+c 



< 



whereas if a € [a/(c + 1), a), 

a a a 

The claim is proven. 

Let Aq G "^[6, cj~^] and define recursively An+i = f{An). We have An £ 
^[6, cr"^] for all n, and for n > 1, as in Section 4, 

.1/2 .-1/2 _ .-1/2 .1/2 _ ^ Al/^^aa*Al/'^a*{An - 



Letting 



we have 



zr _ /.1/2 /.-1/2 _ .-1/2 .1/2 



= cE 



cE 



(1 + a*AnCx){l + aM„_iQ!) 
{l + a*Ana){l + a*An-ia) ' 



or in vector form, 



veci?„+i 



cE 



cE 



(1 + Q*A„Q!)(1 + a*An-ia) 



vecF* 



vecH* 



(l + aM„_i)(l + aM„a) 

using Lemma 2.6. Arguing the same way as in the previous section, we have 
by Lemma 2.8, 



cE 



il + a*An-i){l + a*Ana) 



< 



cE 



Al^^aa*A„aa*Al^^ 



(1 + QM„_ia)(l + a*AnCx) 



cE 



(1 + aM„_ia)(l + qM^q) 
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< 



cE 



1 + a*Ar.-ia 



= - O-^AnW \\Il - O-'^An+lW 
<il-CTHf. 

Notice p=l — a'^b G (0, 1). For n > 2, we therefore get 



l-f^nll < ||veci/„|| < p" ^llvecffill, 



and so 



An, A' 



n-1 1 



Ai/'HnAlZ\\<<T-'p'"'\\yecH,\\, 



which imphes for m>n>2, 



llvecffill 

Therefore, {A^} is a Cauchy sequence, and hence, convergent to a matrix 
A G cj2]. From continuity of /, ^ satisfies (1.1). 



6. Proofs of corollaries. For the first corollary we see that, under as- 
sumption (1.2), / maps diagonal matrices consisting of positive diagonal 
elements into diagonal matrices. Due to the uniqueness of solutions to (1.1), 
the proof amounts to showing the existence of positive ai, . . . ,aL satisfying 
(1.1). This is achieved by invoking Lemma 2.7. We simply take X = M^, g the 
right-hand side of (1.3) and A= [b,a~'^]^. The first statement in Corollary 
1.1 follows. 

For the second corollary, we follow along the argument toward the end 
of Section 3. We see immediately that (3.7) holds. Consider a realization in 
which (3.4), (3.5) and (3.7) hold, the empirical distribution of (|afe(l)p, . . . , 
|Qfc(L)p), 2 < k < K, converges weakly to G, and (1.5) is true for all positive 
rational ai,. . . ,aL. Then, for this realization, a simple continuity argument 
reveals (1.5) true for all positive ai,...,aL- Moreover, the empirical dis- 
tribution oi a^, 2 < k < K , is tight. The subsequence {Ni} considered can 
therefore also be one in which the empirical distribution of ct^, 2 < k < K , 
converges weakly to, say, H. The rest of the argument at the end of Section 
3 leads to only one solution A on {Ni} satisfying (1.1). But, by the domi- 
nated convergence theorem, (1.2) holds for all positive ai, . . . , ai- Thus, from 
Corollary 1.1, A is diagonal satisfying (1.3), which depends only on G. Thus, 
Corollary 1.2 follows. 
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7. Question of uniformity. Let 



N 

Then 



SIRfc = ^PliCk + a^ir'Pk = ^ E ak{i)ak{£')sUCk + a''l)i}sk 

represents user fc's best signal-to-interference ratio. We are interested in 
knowing what conditions are needed to insure 



max 

k<K 



a.s. 



Clearly nothing can be concluded without assuming bounds or some 
growth rate on the afc(^)'s along with knowledge of the rate of convergence 
of the (l/A^)s^(Cfc + (T^/)^|,Sfc's. The latter is tied closely with its limiting 
distributional behavior, which will be investigated in later work. For now we 
will confine the analysis to providing conditions to ensure for any 

(7.1) nMx; |iV~^s^(Cfc + (T^/)^^,Sfc — a^^^'l — > a.s. 

as — > oo. 

We have the following: 



Theorem 7.1. //, in addition to the conditions in Theorem 1.1, E|sii|^ < 
oo, or if the doubly infinite array assumption is dropped, E|sii|^ < oo, then 
(7.1) is true. 

Proof. For each k < K, let s/. denote the vector obtained after trun- 
cating and centralizing s^, the same way as si. Each remains unchanged. 
We have 

\N-^sUCk + a^I)i}sk - N-%{Ck + <t2/)-i sfcl 

< a~^N~W\sk - Skf + 2||sfc||||sfc - Sfell). 

By Lemma 2.11, {Xu = |siip, a = /3 = 1) under the double array assump- 
tion, or, for nondouble array. Lemma 2.10 {Xi = t = 3) together with 
Boole's inequality, we follow the steps in the beginning of Section 3 and find, 
almost surely, 

limsupmax I A^~^ l|si.|P — 1| = 
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and 

limsupmax A^~^ ||sfc — s^lp = 0. 

N k<K 

Letting = Sfc/(E|siip)"'^/^, it follows that, almost surely, 

maxiV-^|s^(C7fc + a'^I)^},Sk - sl{Ck + a'^I)llsk\ ^ 

as — > oo. Applying Lemma 2.9 for any p > 4, we have then 

maxN^^\sl{Ck + (T^/)7isfc — tr(Cfc + a'^I)j], \ a.s. 

k<K ' ' 

as N ^ oo. For any k,k' , with two applications of Lemma 2.1, we find (with 

Thus, the remainder of the proof of Theorem 1.1 proceeds exactly as in 
Section 3, and we get (7.1). □ 
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